5.3 Methods for Analysis and Processing of Discrete Biosignals

|

199

about half. After that, the amplification increases again, since the spectrum repeats

according to the sampling theorem (cf. Figure 5.5). This effect is due to the fact that

the analogue RC low-pass filter used does not have an infinitely large attenuation at

the cut-off frequency fg like an ideal low-pass filter, but only an attenuation of 3 dB,

which is equivalent to a drop in the magnitude square of the transfer function by

1/20.707. At higher frequencies, the attenuation increases but is non-zero at half

the sampling frequency of fab/2 = 500 Hz. This means that there is overlapping of the

spectrum after sampling. This cannot result in an ideal low-pass filter that has an in-

finitely large attenuation at its cut-off frequency. Furthermore, the property of the cor-

responding analogue filter according to Figure 5.39, which also has an infinitely large

attenuation at a higher frequency. A remedy is achieved by not requiring the same val-

ues as those of the impulse response of the analogue filter for the sampling times in

the filter design, but only striving for the closest possible approximation of the trans-

fer function. However, because of the basically periodic course of the transfer function

of the digital system, this means that the transfer function of the analogue system in

the range from 0 tocorresponds to that of the digital system in the range from 0

to half the sampling frequency. The frequency axis of the analogue system is virtually

mapped onto the frequency axis of the digital system from 0 to half the sampling fre-

quency. This can be done according to the figure after the bilinear transformation with

the following procedure:

1.

The transfer function G(z = ejωTa) shall match that of the associated analogue

filter G(p = up to half the sampling frequency, i.e..

G(z = ejωTa) = G(p =) .

(5.93)

2.

substituting the inverse function of z, i.e..

p = ln(z)

Ta

,

(5.94)

into the transfer function G(p) of the analogue filter does not lead to a fractional

rational function in z, and the digital filter is thereby not realisable with adders,

delayers and multipliers (as, for example, in the canonical circuits).

3.

The ideal transformation between p- and z-range, z = epTa or p = ln(z)

Ta , is therefore

approximated in such a way that a realisable transfer function in the z- range is

obtained.

4.

To do this, ln(z) is generated by a series expansion with fractional functions in

z and broken off after the 1st member. This is then an approximation to be used

called bilinear transformation; bilinear here means that there are linear functions

in p or z in both the numerator and denominator:

ln z = 2 { z1

z + 1 + 1

3 ( z1

z + 1)

3

+ 1

5 ( z1

z + 1)

5

+ ⋅⋅⋅}2 z1

z + 1 .

(5.95)